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^ ■ Abstract. We study the Polyakov loop dynamics originating from finite-temperature Yang- 

Mills theory. The effective actions contain center-symmetric terms involving powers of the 
Polyakov loop, each with its own coupling. For a subclass with two couplings we perform 
iy-j | a detailed analysis of the statistical mechanics involved. To this end we employ a modified 

mean field approximation and Monte Carlo simulations based on a novel cluster algorithm. 
| We find excellent agreement of both approaches. The phase diagram exhibits both first and 

second order transitions between symmetric, ferromagnetic and antiferromagnetic phases 
with phase boundaries merging at three tricritical points. The critical exponents v and 7 at 
■ the continuous transition between symmetric and antiferromagnetic phases are the same as 

C f° r the 3-state spin Potts model. 
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' tion; Monte Carlo simulations 
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■ 1 Introduction 



Symmetry constraints and strong coupling expansion for the effective action describing the 
Polyakov loop dynamics of gauge theories lead to effective field theories with rich phase struc- 
tures. The fields are the fundamental characters of the gauge group with the fundamental 
domain as target space. The center symmetry of pure gauge theory remains a symmetry of the 
effective models. If one further freezes the Polyakov loop to the center Z of the gauge group 
one obtains the well known vector Potts spin-models, sometimes called clock models. Hence 
we call the effective theories for the Polyakov loop dynamics generalized Z-Potts models. We 
review our recent results on generalized Z3-Potts models pQ. These results were obtained with 
the help of an improved mean field approximation and Monte Carlo simulations. The mean field 
approximation turns out to be much better than expected. Probably this is due to the existence 
of tricritical points in the effective theories. There exist four distinct phases and transitions of 
first and second order. The critical exponents v and 7 at the second order transition from the 
symmetric to antiferromagnetic phase for the generalized Potts model are the same as for the 
corresponding Potts spin model. 

Earlier on it had been conjectured that the effective Polyakov loop dynamics for finite 
temperature SU(N) gauge theories near the phase transition point is very well modelled by 
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3-dimensional Z^r spin systems OS]- For SU(2) this conjecture is supported by universality 
arguments and numerical simulations. The status of the conjecture for SU(3) gauge theories is 
unclear, since the phase transition is first order such that universality arguments apparently are 
not applicable. 



2 Recall of planar Potts models 

The g-state Potts model [UEJE] is a natural extension of the Ising model. On every lattice site x 
there is a planar vector with unit length which may point in q different directions (Fig. [1]) with 
angles 




2tt 4-7T 
q q 



,2ir 



2, 



Only nearest neighbors interact and their contribu- 
tion to the energy is proportional to the scalar prod- 
uct of the vectors, such that 



H 



- J cos 



(2.1) 



The Hamiltonian H is invariant under simultaneous 
rotations of all vectors by a multiple of 2k jq. These 
Figure 1. Angles in a g-state Potts model. Z « symmetries map a configuration w = {9 x \x € A} 

into w' = {8 X + 2irn/q}, n £ {1, . . . , q}. In two and 
higher dimensions the spin model shows a phase transition at a critical coupling K c = (3J C > 
from the symmetric to the ferromagnetic phase. In two dimensions this transition is second 
order for q < 4 and first order for q > 4. In three dimensions it is second order for q < 2 and 
first order for q > 2. 

There is another phase transition at negative cri- 
tical coupling K' c from the symmetric to the anti- 
ferromagnetic phase. For q > 3 and negative K the 
number of ground states of (|2.ip increases rapidly 
with the number of lattice sites. It equals the num- 
ber of ways one can color the vertices with q colors 
such that two neighboring sites have different colors 
(Fig. [2]). In the antiferromagnetic case the dege- 
nerate ground states contribute considerably to the 
entropy [7] 




Figure 2. Coloring of neighboring lattice 
sites. 



s B (p) = -J2 p ( w ) l °g p ( w )> 



where one sums over all spin configurations and P is 
the probability of w. We use the well known varia- 



tional characterization of the free energy, 

(3F = inf (0(H) P - S B ) , (H) P = ]T P(w)H(w), 

w 

where the minimum is to be taken on the space of all probability measures. The unique mini- 
mizing probability measure is the Gibbs state 



Pi 



Gibbs ~ e 



-f3H 



(2.2) 



Generalized Potts-Models and their Relevance for Gauge Theories 



3 



belonging to the canonical ensemble. In the variational definition of the convex effective action 
one minimizes on the convex subspace of probability measures with fixed mean field, 

r[m] =mf(P(H)p-S(P)\(e ie *) P = m(x)). (2.3) 

The field m(x) which minimizes the effective action is by construction the expectation value of 
the field e x in the thermodynamic equilibrium state (j2.2|) . 

In the mean field approximation to the effective action one further assumes that the measure 
is a product measure [HIE], 

P(w) = P({8 x }) = l[p x (9 x ), (2.4) 

X 

where the single site probabilities are maps p x : Z q — > [0, 1] . The approximate effective action is 
denoted by r m f[m]. 

The symmetric and ferromagnetic phases are both translationally invariant. In the mean field 
approximation all single site probabilities are the same, p x = p, and the mean field is constant, 
m(x) = m. The effective potential is the effective action for constant mean field, divided by the 
number of sites. Its mean field approximation is 

u mf (m) = inf {-Kmm* + ^p{6)\ogp{6)\^p{6)e ie = m\ , K = dJ. (2.5) 
P V 6 e J 

It agrees with the mean field approximation to the constraint effective potential, introduced by 
O'Raifeartaigh et al. [101 111] • In the antiferromagnetic phase there is no translational invariance 
on the whole lattice A but on each of two sublattices in the decomposition A = Ai U A2. The 
sublattices are such that two nearest neighbors always belong to different sublattices. Thus the 
single site distributions p x in (|2.4h are not equal on the whole lattice, but only on the sublattices, 

Px = Pi on Ai and p x = P2 on A2. (2-6) 

The minimization of the effective action on such product states is subject to the constraints 

pi(6)e ie = mi and P2{0)e ie = m 2 

0G Zq 06 Zq 

and yields the following mean field effective potential 

u mf (mi, m 2 ) = i ^K\mi - m 2 | 2 + ^ u mf (mj)j , 

where u m { is the effective potential (|2.5|) . For K > the minimum is attained for mi = m-2 
and translational invariance is restored. In the symmetric and ferromagnetic phases the single 
site probabilities p\ = p 2 = P- In the symmetric phase p = 1/q for every orientation and in the 
ferromagnetic phase p is peaked at one orientation. Hence there exist q different ferromagnetic 
equilibrium states related by Z q symmetry transformations. In the antiferromagnetic phase the 
probabilities p\ and p 2 are different. 

For the 3-state Potts model one can calculate the single site probabilities explicitly. On one 
sublattice it is peaked at one orientation and on the other sublattice it is equally distributed 
over the remaining 2 orientations. Hence there are 6 different antiferromagnetic equilibrium 
states related by Z3 symmetries and an exchange of the sublattices. The results for single site 
distributions of the 3-state Potts model are depicted in Fig. [3j 
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A 





symmetric ferromagnetic antiferromagnetic 

Figure 3. Single site distributions of the 3-state Potts model. 



3 Polyakov-loop dynamics 

We consider pure Euclidean gauge theories with group valued link variables U x -n on a lattice with 



Nt sites in the temporal direction. The fields are periodic in this direction, U 



t+N-, 



u 



We are interested in the distribution and expectation values of the traced Polyakov loop variable 



Ni 



L x = trV x , V x = JJ U ttX - { 



(3.1) 



t=i 



since (L x ) is an order parameter for finite temperature gluodynamics (see |12j for a review). 
In the low-temperature confined phase (L x ) = and in the high-temperature deconfined phase 
(L x ) 0. The effective action for the Polyakov loop dynamics is 



-SenlV] 



VU5 




(3.2) 



links 



with gauge field action S w . In this formula the group valued field V x is prescribed and the delta- 
distribution enforces the constraints (|3.ip . In the simulations we used the Wilson action for the 
gauge fields. Gauge invariance of the action S W [U] and measure T>U implies SWI^ 7 ] = 5 e ff[L]. 
In addition there is the global center symmetry, under which all V x are multiplied by the same 

center element of the group. For SU (N) the center 
consists of the N roots of unity, multiplied by the 
identity matrix. Hence for SU(N) theories we have 




S eS [L] = S cf[ [z ■ L], 



1. 



Figure 4. Fundamental domain of L. 



In Fig. [J] on the left we plotted the domain of the 
traced Polyakov loop variable for SU (3). The values 
of L at the three center elements are 3, 3z and 3z 2 
with z = e 27 ™/ 3 . They form the edges of the trian- 
gle. What is needed is a good ansatz for the effective 
action S e g in f|3.2[) . To this aim we calculated the 
leading terms in the strong coupling expansion for 
S e g in gluodynamics [TJ [13] . As expected one finds 
a character expansion with nearest neighbor inter- 
actions, 

3N t \ 



Seff — AioSio + A20S20 + AnSii + A21S21 + P1V1 + 0((3 

where S pq depends on the character Xpg belonging to the representation (p, q) of 577(3). Ex- 
pressing the characters as function of the fundamental characters xio = L and xoi = L* the 
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different center-symmetric contributions have the form 



Sio 


= V (L X L* + h.c.) , 
(xy) 








= V ( L 2 L* 2 - L 2 L„ - 

/ j y-'x^y ^x^y 
(xy) 


x y 1 x y 


+ h.c.) 




= ^] ( \L X 2 \Ly 2 — ij. 

(xy) 


2 -\L y \ 2 + l), 




S21 




2L* x L y + h.c.) , 


V x 



= E(i^i 2 - 1 )- 



The target space for L is the fundamental domain inside the triangle depicted above. The 
functional measure is not the product of Lebesgue measures but the product of reduced Haar 
measures on the lattice sites. 

In the following we consider the two-coupling model 

Seff = (Aio - 2A 21 ) (L x L* y + h.c.) + A 21 (L 2 x L y + L 2 y L x + h.c.) (3.3) 

which contains the leading order contribution in the strong coupling expansion. For vanishing A21 
it reduces to the Polonyi-Szlachanyi model [T7] . 



4 Gluodynamics and Potts-model 

There is a direct relation between the effective Polyakov loop dynamics and the 3-state Potts 
model: If we freeze the Polyakov loops to the center of the group, 

f 2ir 47r 1 

V x — ► ZjlG center (SU(3)) 8 X G i 0, — , — \ 

then the effective action (13. 3p reduces to the Potts-Hamiltonian (|2.ip with J = 18(Aio + 4A21). 
The same reduction happens for all center-symmetric effective actions with nearest neighbor 
interactions. Only the relation between the couplings X pq and J is modified. 

For the gauge group SU(2) the finite temperature phase transition is second order and the 
critical exponents agree with those of the 2-state Potts spin model which is just the ubiquitous 
Ising model. The following numbers are due to Engels et al. [14j 





(3/u 7/1/ v 


Ad SU(2) 
3d Ising 


0.525 1.944 0.630 
0.518 1.970 0.629 



and support the celebrated Svetitsky-Yaffe conjecture [2_p. For the gauge group SU(3) the finite 
temperature phase transition is first order and we cannot compare critical exponents. But the 
effective theory (|3.3p shows a second order transition from the symmetric to a antiferromagnetic 
phase and we can compare critical exponents at this transition with those of the same transition 
in the 3-state Potts spin model. 

In a first step we study the 'classical phases' of the Polyakov loop model with action (|3.3p . 
The classical analysis, where one minimizes the 'classical action' S^g, shows a ferromagnetic, 
antiferromagnetic and anticenter phase. The classical phase diagram is depicted in Fig. For 



For the relations between 3-dimensional gauge theories at the deconfining point and 2-dimensional Potts- 
models, the so-called gauge-CFT correspondence, we refer to the recent paper [15] . 
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Figure 5. Classical phase diagram. 



small couplings the quantum fluctuations will disorder the system and entropy will dominate 
energy. Thus we expect a symmetric phase near the origin in the (Aio, A2i)-plane and such 
a phase was inserted by hand in the diagram. There exists one unexpected phase which cannot 
exist for Potts spin models. It is an ordered phase for which the order parameter is near to the 
points opposite to the center elements, hence we call it anticenter phase. They are marked in the 
fundamental domain in Fig. 21 We shall see that the classical analysis yields the qualitatively 
correct phase diagram. 

4.1 Modified mean field approximation 

In a second step we calculate the effective potential for the Polyakov loop model with action (|3.3p 
in the mean field approximation. Here we are not concerned with the relevance of this model 
for finite temperature gluodynamics and just consider (j3 . 3|) as the classical action of a field 
theoretical extension of the Potts spin models. We use the variational characterization (|2.3p 
for the effective action where we must minimize with respect to probability measures on the 
space of field configuration {V x \x G A} with fixed expectation values (xpq) of all characters 
Xpqi'Px) showing up in the Polyakov loop action. As outlined in Section [2j in the mean field 
approximation we assume the measures to have product form, 



For further details the reader is referred to our earlier paper |16j . Here /i re( j is the reduced Haar 
measure of SU(3). Since we expect an antiferromagnetic phase we only assume translational 
invariance on the sublattices in the decomposition (|2.6p . This way one arrives at a non-trivial 
variational problem on two sites. 

We illustrate the procedure with the simple Polyakov loop model studied in |17j . 



To enforce the two constraints (L x ) = Li for x G A, one introduces two Lagrangian multipliers. 
For the minimal model (|4.ip one arrives at the following mean field effective potential 




(4.1) 



2u m i(L 1 ,L* 1 ,L 2 ,L* 2 ) = -dX^ - L 2 \ 2 + Y,Vmi(Li,L*) 
with v m{ (L,L*) = d\\L\ 2 + j (L,L*). 
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(a) 1.6 



1.4 
1.2 
1.0 ■ 

I 

' 0.8 ■ 
0.6 
0.4 
0.2 ■ 



Aio, c = -0.13433(1) 
At = 1.46{2) 



mean field 



o.o I 



-0.1350 



-0.1345 



-0.1340 



(b) 



1.6 
1.4 

1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
o.u 



Aio >c - -0.13721(5) 
M ~ 1-33(2) 



multicanonkal 



Monte Carlo 



-0.1335 



0.25 



0.20 



|0.15 

m 

0.10 



0.05 



-0.1378 



-0.1374 



-0.1370 



Figure 6. The order parameter i for (a) mean field approximation and (b) Monte Carlo simulation. 
Here 70 is the Legendre transform of 

Wo(j,f) = log0 O (i,i*), z (j,f) = J dfi rcd exp (jL + j*L*) . 

The last integral has an expansion in terms of modified Bessel functions [181 119j . 



inNa 



In 


In+1 


In+2\ 






In-l 
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In+1 (2|i|) , 


J = \ 


j\e ia 


In-2 


In-l 


In J 







As order parameters discriminating the symmetric, ferromagnetic and antiferromagnetic phases 
we take |20] 



L = ±(L 1 + L 2 ) 



M=^(L 1 -L 2 ), 



\L\ 



m = \M\ 



(4.2) 



Fig. [6] shows the value of the order parameter t as function of the coupling Ain near the phase 
transition point from the symmetric to the ferromagnetic phase. The value of the critical cou- 
pling and the jump A£ of the order parameter in the mean field approximation and simulations 
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0.188 



Figure 7. The order parameter m for (a) mean field approximation and (b) Monte Carlo simulation. 



agree astonishingly well. In Fig. [6j3 the order parameter £ for the ferromagnetic phase is plotted 
against its probability distribution given by the shaded area. In the vicinity of the first order 
transition we used a multicanonical algorithm on a 16 3 lattice to calculate the critical coupling 
to very high precision. For further details on algorithmic aspects I refer to our recent paper pQ. 
Since the model has the same symmetries and dynamical degrees as the 3-state Potts model 
in 3 dimensions we expect the model to have a second order transition from the symmetric 
to the antiferromagnetic phase. To study this transition we calculated the order parameter m 
?2j) in the modified mean field approximation and with Monte Carlo simulations. The order 



m 



parameter m is sensitive to the transition in question. 

Again the expectation value and probability distribution of m near the transition is plotted 
in Fig. [7J To get a clear signal we have chosen a large lattice with 28 3 sites and evaluated 
5 x 10 5 sweeps. The mean field approximation and Monte Carlo simulations together with 
the calculation of Binder cumulants in [1] indicate that the transition to the antiferromagnetic 
phase is second order. We plan to use a finite size scaling method to confirm this result in our 
upcoming work. 

With the cumulant method we have calculated the critical exponents 7 and v and compared 
our results with the same exponents for the 3-state Potts model at the second order transition 
from the symmetric to the antiferromagnetic phase [20] . Within error bars the critical exponents 
are the same. 
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exponent 


3-state Potts minimal S e s 


V 

7f/v 


0.664(4) 0.68(2) 
1.973(9) 1.96(2) 



This is how the conjecture relating finite temperature gluodynamics with spin models is at 
work: The Polyakov loop dynamics is effectively described by generalized Potts models with the 
fundamental domain as target space, for SU{3) it is the triangularly shaped region in Fig. 01 The 
first order transition of gluodynamics is modelled by the transition from the symmetric to the 
ferromagnetic phase in the generalized Potts models. These generalized Potts models are in the 
same universality class as the ordinary Potts spin models; they have the same critical exponents 
at the 'unphysical' second order transition from the symmetric to the antiferromagnetic phase. 

It is astonishing how good the mean field approximation is. The reason is probably, that the 
upper critical dimension of the (generalized) 3-state Potts model is 3 and not 4 as one might 
expect. This is explained by the fact, that the models are embedded in systems with tricritical 
points, see below, and for such systems the upper critical dimension is reduced |21^ I22j. 



5 Simulating the effective theories 

We have undertaken an extensive and expensive scan to calculate histograms in the coupling 
constant plane (Aio, A21) of the model (13. 3p . Away from the transition lines we used a standard 
Metropolis algorithm giving results within 5 percent accuracy. Near first order transitions 
we simulated with a multicanonical algorithm on lattices with up to 20 3 lattice sites. Most 
demanding have been the simulations near second order transitions. For that we developed a new 
cluster algorithm [T] which improved the auto-correlation times by two orders of magnitude on 
larger lattices. We found a rich phase structure with 4 different phases with second and first 
order transitions and tricritical points. As for the minimal model the Monte Carlo simulations 
are in good and sometimes very good agreement with the mean field analysis. 

Fig. [8] shows the phase structure in the generalized MF approximations and the corresponding 
results of our extended MC simulations. The results of the simulations are summarized in the 
following phase portrait (Fig. [9]), in which we have indicated the order of the various transitions. 
The calculations were done on our Linux cluster with the powerful jenLaTT package. In 3000 
CPU hours we calculated 8000 histograms in coupling constant space. 




-.2 -.1 0.1 0,2 0.3 -.2 -.1 0.1 0.2 0.3 



Figure 8. Phase structure for MF approximation and MC simulation. 
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Aio = -0.13958 
A 2 i = 0.0020833 



= -0.13971 
■ 0.0019583 




Aio = -0.13996 
A 2 i = 0.0017083 





Aio = -0.13983 
A21 = 0.0018333 



A I0 = -0.14008 
A21 = 0.0015833 





Aio = -O.14021 
A21 = 0.0014583 




Figure 10. Histograms of L which show a clear signal for the first order transition from the symmetric 
to the ferromagnetic phase. 

For each marked transition in Fig. [9] a set of 6 histograms is displayed. These histograms 
(and many others, see pQ) have been used to localize the phase transition lines and to investigate 
the nature of the transitions. The results are summarized in the phase portrait on page 1101 The 
critical exponents v and 7 given above have been determined at the second order transition 
indicated with arrow [13] in the portrait. 
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Figure 11. Histograms of L for the continuous transition from the symmetric to the ferromagnetic 
phase. 
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Figure 12. Histograms of L for the continuous transition from the symmetric to the anticenter phase. 
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Figure 13. Histograms of M for the second order transition from the symmetric to the antiferromagnetic 
phase. 



6 Conclusion 



The strong coupling expansion results in a character expansion for the Polyakov-loop dynamics. 
The leading terms are center symmetric nearest neighbor interactions containing the characters 
of the smallest representations of the gauge group. We have performed an extensive modified 
mean field analysis which includes antiferromagnetic states without translational invariance on 
the whole lattice. A new and efficient cluster algorithm has been developed and applied to 
study the second order transitions from the symmetric to the antiferromagnetic phase. The 
autocorrelation times were improved by 2 orders of magnitude. We discovered an unexpectedly 
rich phase structure of the simple 2-coupling Polyakov loop model (|3.3p . This model is in the 
same universality class as the 3-state Potts spin model. The mean field results are surprisingly 
accurate that seems to indicate that the upper critical dimension of the generalized Potts models 
is 3. This is attributed to the existence of tricritical points [2T] [22] . 

To relate our results to gluodynamics we have to calculate the effective couplings X pq gover- 
ning Polyakov-loop dynamics as functions of the Wilson coupling (3 in gluodynamics. We have 
done this successfully for SU (2) gauge theory with inverse Monte Carlo techniques |16[ [23] and 
plan to publish our results for SU(3) very soon [24] . For the inverse Monte Carlo simulations to 
work one needs simple geometric Schwinger Dyson equations for the Polyakov loop dynamics. 
Such equations have been derived very recently in [19]. It would be interesting to see whether 
the antiferromagnetic phase of the Polyakov loop models plays any role at all for gluodynamics. 
In the present paper it was needed to show that certain critical exponents of the Polyakov loop 
models are the same as of the q = 3 Potts spin model. Finally one would like to include heavy 
fermions in the effective Polyakov-loop dynamics [25]. To that end one needs to add center 
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symmetry breaking terms to the effective actions studied in the present paper. This will lead 
to a proliferation of additional terms in the effective action which renders a systematic study 
more difficult as compared to pure gluodynamics. 
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